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Abstract 

A review of a new separability theory based on degenerated Poisson pencils and the so- 
called separation curves is presented. This theory can be considered as an alternative 
to the Sklyanin theory based on Lax representations and the so-called spectral curves. 

1 Introduction 

The separation of variables belongs to the basic methods of mathematical physics from 
the previous century. Originating from the early works of D'Alambert (the ISth century), 
Fourier and Jacobi (the first half of the 19th century), for many decades it has been 
the only known method of exact solution of dynamical systems. Let us briefly recall the 
idea from classical mechanics. Consider a Hamiltonian mechanical system of 2n degrees 
of freedom and integrable in the sense of Liouville/Arnold theorem. This means that 
there exist n linearly independent functions hi (Hamiltonians) which are in involution with 
respect to the canonical Poisson bracket 

{/ij,/ifc} = 0, j,k = l,...,n. (1.1) 
A system of canonical variables (//j, Xi)^=i 

{Xi, Xj} = {fXi, fij} = 0, {Xi,^ij} = 6ij (1.2) 
is called separated if there exist n relations of the form 

<^i{Xi,tii,hi,...,hn) =0, i = l,...,n (1.3) 

joining each pair (/ij,Aj) of conjugate coordinates and all Hamiltonians h^, k = l,...,n. 
Fixing the values of Hamiltonians hj = const j = aj one obtains from (|1.3| ) an explicit 
factorization of the Liouville tori given by the equations 

<^i{Xi,Hi,ai,...,an) = 0, i = l,...,n. (1.4) 
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In the Hamilton- Jacobi method, for a given Hamiltonian function hr{n, A) we are looking 
for a canonical transformation (/U, A) {(i,b) in the form bi = ^-^IJ-i = fy-, where 
W{\,a) is a generating function given by the related Hamilton-Jacobi equation 

dW 

hr[X,—) = ar. (1.5) 
If (/i, A) are separated coordinates, then 

n 

W{\a) = Y,W,{\ua) (1.6) 

i=l 

and the partial differential equation (^]^) splits into n ordinary differential equations 
dW 

¥'j(Ai, -^,ai, ...,an) = 0, ? = l,...,n (1.7) 



just of the form ( |1.4[ ). In (a, b) coordinates the flow is trivial 

{a,)u = 0, {bX = Sjr (1.8) 
and the implicit form of the trajectories Xi{tr) is the following 
dW 

bj(X,a) = — — = Sjrtr + const, j = l,...,n. (1-9) 

oaj 

So, given a set of separated variables, it is possible to solve a related dynamical system by 
quadratures. In the 19th century and most of the present century, for a number of models 
of classical mechanics the separated variables were either guessed or found by some ad hoc 
methods. For example, in the second half of the 19th century, Neumann's investigation of 
a particle moving on a sphere under the action of a linear force |Q] and Jacobi's study of 
the geodesies motion on an ellipsoid exploited the separability of the Hamilton-Jacobi 
equation to solve the equations of motion by quadratures. In 1891 Stackel initiated the 
program dealing the classification of Hamiltonian systems according to their separability 
or nonseparability, presenting conditions for separability of the Hamilton-Jacobi equation 
in orthogonal coordinates Q. For three dimensional flat space, the 11 possible coordinate 
systems in which separation may take place were deduced in a paper by Eisenhart [^. 
They were all obtained as degenerations of the confocal ellipsoidal coordinates ||6|. For 
each of the coordinates Eisenhart Q determined the form of the potential that permitted 
a separation of variables. These potentials, designated Stackel or separable potentials, 
played a crucial role in Hamiltonian mechanics before the development of more qualitative 
geometric methods for differential equations. Although in all classical papers the trans- 
formation to separated coordinates was searched in the form of point transformations, 
nevertheless they can be produced by an arbitrary canonical transformation involving 
both coordinates and momenta. 

A fundamental progress in the theory of separability was made in 1985, when Sklyanin 
adopted the method of soliton systems, i.e. the Lax representation, to systematic deriva- 
tion of separated variables (see his review article [Q] ) . It was the first constructive theory 
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of separated coordinates for dynamical systems. In his approach, the appropriate Hamil- 
tonians appear as coefficients of the spectral curve, i.e. the characteristic equation of Lax 
matrix. His method was successfully applied to separation of variables for many old and 
new integrable systems [p[|- [p^] . 

In this paper we present a new constructive separability theory, which has been re- 
cently intensely developed, based on a bi-Hamiltonian property of integrable systems. In 
last decade a considerable progress has been made in construction of new integrable fi- 
nite dimensional dynamical systems showing bi-Hamiltonian property. The majority of 
them originate from stationary flows, restricted flows or nonlinearization of Lax equations 
of underlying soliton systems |13|-p5|. Quite recently a fundamental property of such 
systems has been discovered, i.e. their separability. It was proved |^]-|31] that most 
bi-Hamiltonian finite dimensional chains, which start with a Casimir of the first Poisson 
structure and terminate with a Casimir of the second Poisson structure, are integrable by 
quadratures, through the solutions of the appropriate Hamilton-Jacobi equation. 

The presented review article is based on results of papers p6|-[^] systematizing and 
unifying them into a compact separability theory in the frame of the set of canonical 
coordinates. The results derived so far are sufficiently promising to consider the theory as 
an alternative or complement to the Sklyanin ones. 



2 Preliminaries 

Let us re-examine some facts about bi-Hamiltonian systems. We recall some definitions. 
Let M be a differentiable manifold, TM and T*M its tangent and cotangent bundle. 
At any point u E M, the tangent and cotangent spaces are denoted by T^M and T*M, 
respectively. The pairing between them is given by the map < •, • >: T*M x T„M R. 
For each smooth function F £ C°°{M), dF denotes the differential of F. M is said to be 
a Poisson manifold if it is endowed with a Poisson bracket {•, •} : C°°(M)x C°°{M) — > 
C°°(M), in general degenerate. The related Poisson tensor vr is defined by {F, G}t^{u) :=< 
dG,TT o dF > (u) =< dG{u),'ir{u)dF{u) >. So, at each point u, 7r(u) is a linear map 
7r{u) : T*M — > T^M which is skew-symmetric and fulfils the Jacobi identity. Any function 
c G C°°(M), such that dc £ kervr, is called a Casimir of vr. Let 7ro,vri : T*M TM be 
two Poisson tensors on M. A vector field K is said to be a bi-Hamiltonian with respect to 
ttq and vri if there exist two smooth functions H,F € C°°(M) such that 

K = 7rQodH = 7rio dF. (2.1) 

Poisson tensors ttq and vri are said to be compatible if the associated pencil ttx = tti — Xttq is 
itself a Poisson tensor for any X £ R. Moreover, if ttq is invertible, the tensor = vri ottq ^, 
called a recursion operator, is a Nijenhuis (hereditary) tensor of such a property that when 
it acts on a given bi-Hamiltonian vector field K, it produces another bi-Hamiltonian vector 
field being a symmetry generator of K. Hence, having the invariant Nijenhuis tensor, one 
can construct a hierarchy of Hamiltonian symmetries and related hierarchy of constants 
of motion for an underlying system, so important for its integrability. 

Unfortunately, for majority of bi-Hamiltonian finite dimensional systems, both Poisson 
structures are degenerate, so one cannot construct the recursion Nijenhuis tensor inverting 
one of the Poisson structures. Nevertheless, due to the nonuniqueness of Hamiltonian 
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functions, determined up to an appropriate Casimir function, it is always possible to 
construct a finite bi-Hamiltonian chain starting and terminating with Casimirs of ttq and 
vTi, respectively. 

3 One- Casimir chains 

Let us consider a Poisson manifold M of dim M = 2n + 1 equipped with a linear Poisson 
pencil 

TTA = TTl - AtTo (3.1) 

of maximal rank, where ttq and vri are compatible Poisson structures and A is a continuous 



parameter. As was first shown by Gel'fand and Zakharevich |32|, a Casimir of the pencil 
is a polynomial in A of an order n 

/iA = /ioA" + /iiA"-i + ... + /i„ (3.2) 

and generates a bi-Hamiltonian chain 

ttq o dhQ = 

ttq o dhi = Ki = TTl o dhQ 

TTo O dh2 = K2 = TTl o dhi 

TTx o dhx = <^=^ . (3-3) 

TTq O dhn = Kn = TTl o dhn-l 
= TTl o dhni 

where K = Ki, H = /iiand F = hQ.li). the paper we restrict our considerations to a class 
of canonical coordinates {q,p,c), where q = {qi, ...,qn)'^ , p = {pi, ...,Pn)'^ are generalized 
coordinates and c is a Casimir coordinate. Hence, vro always stays a canonical Poisson 
matrix. This restriction simplifies the theory in the sense that it makes a Marsden-Ratiu 
projection procedure |3^], |Q trivial. 



3.1 Darboux-Nijenhuis representation 

To understand the theory better, let us start from the end in some sense, i.e. from a 
system written in separated coordinates (^j,Aj)f^]^. An interesting observation is that 
such a system can be represented by n different points of some curve. We start from 
Gel'fand and Zakharevich case. Actually, let us consider a curve in (A, ^) plane, in the 
particular form 

fix, fi) = hx, hx = cA" + /iiA"-i + ... + K, (3.4) 

where /(A,/x) is an arbitrary smooth function. Then, let us take n different points (/Xj, Aj) 
from the curve: 

/(A„/xJ =cA^ + /ii A^-i + ... + /!„, i = l,...,n, (3.5) 



which will define our separated coordinates as they are of the form (1.3). The explicit 



dependence of on (//j, Aj, c)f^i is given by the solution of n linear equations ( p. 51) , while 
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for fixed values of hk = ak and fi^ = the system ( |3.5D allows us to solve the appropriate 
Hamilton- Jacobi equations. 

Here, we have to mention that the idea to relate the multi-Hamiltonian property to an 
m-parameter family of curves comes from P. Vanhaecke . 

In refs. p^ ]-p9[| the bi-Hamiltonian chain was constructed for a separation curve in 
the form (|3.4|) . Actually, the Hamiltonian functions hk found from the system (|3.5|) take 
the following compact form 

/ifc(A,^,c) = rvrf- + cpfc(A), k = l,...,n 

i=i ^i^^J (3-6) 

ho = c, 

where 

MX):=ll{X,-Xj), (3.7) 

p,{X) := i-l)' ^ A,,-. ..-A,,, k = l,...,n, (3.8) 

ji, ■■,jk 
ji < ■• < jfc 

are the so-called Viete polynomials (symmetric polynomials) and 

pUW ■■= Pk-iih = 0) = (3.9) 
For example, for n = 2 we have 

Pi = -Xi - A2, P2 = A1A2 (3.10) 
and for n = 3 

Pi = -Xi - A2 - A3, P2 = A1A2 + A1A3 + A2A3, P3 = -A1A2A3. (3.11) 



The bi-Hamiltonian chain ( |3.3| ) is constructed with respect to the following compatible 
Poisson matrix 

A 
-A 

{hi,^)^ (/ii,a) 





( 5 


I 

























) 






A -/ii,A , (3.12) 



where A = diag{Xi, A„) and /ii^^ := (^f^, f^j • Notice that the last column of vri 
is just the first vector field Ki. All Hamiltonians hk are in involution with respect to both 
Poisson structures ttq and vri. 

Applying the following important relations 

\n 

PrW = -Y.Pr-lf^, (3.13) 

4 = 1 

^ dpr Ar _ A P;-i(A) Ar _ [ 1, rn = n-r\ 

Z^dX.\-^ A, -\0, m/n-r /'^-''••'" ^^'^^^ 
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and the decomposition (1^), the Hamilton- Jacobi equations 



dW 

hr{X,-^,c) = ar, r = l,...,n (3.15) 

turn into the form 

^ pti(A)[/(Afc,0VFfc/9Afc) - cXl] _ _ 

/ J ^ — CLr, r— i,...,n (^o.ioj 

k=l ^ 

with the solution 

f{\k,dWk/dXk,c)=cXl + aiXl-^ + ... + an, k = l,...,n. (3.17) 

Hence, W{X, a) can be obtained by solving n decoupled first-order ODEs ( |3.17| ) and the 
family of dynamical systems ( |3.3D can be solved by quadratures. 

Now, let us pass to the projection of the Poisson pencil ttx onto a symplectic leaf S 
of TTo (dim S = 2n) fixing the value of c. Generally, one has to apply the Marsden-Ratiu 
theorem which in this case is trivial, as obviously 6\ = 6i — XOq, where 

is a nondegenerate Poisson pencil on S. Hence, the related Nijenhuis tensor 

iV = eio0-i=(^ °) (3.19) 

is diagonal and this is the reason why we will refer below to the separated coordinates as to 
the Darboux-Nijenhuis (DN) coordinates. Notice that ( p.8| ) are coefficients of minimal 
polynomial of the Nijenhuis tensor 

n n 

{det{N - A))i/2 = A" + ^ ftA"-^ = fJ(A - A^). (3.20) 

i=l i=l 



S the chain (|3.3[ 


turns into the form 






6*0 o dhi = Ki 


= 9i dhn 


Pn 


9o o dh2 = K2 


= ——61 dhn + Oi° dhi 


Pn 


9o dhn = Kn 


= — ^-^^—^61 dhn + OlO dhn 




Pn 



(3.21) 



The last equation terminates the sequence of vector fields Kj. in the hierarchy as for the 
next equation from the chain we have 

6*0 o dhn+i = Kn+i = —Oi o dhn + Gi ° dhn = 0. (3.22) 
Obviously is not a recursion operator for the hierarchy ( ^.21 ). Because of the form of 



the first equation, the vector field Ki is called a quasi-bi-Hamiltonian |36|-[^] and the 



chain ( p. 21 ) could be treated alternatively as a starting point of the separability theory 



for the case of c = 0. 
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3.2 Arbitrary canonical representation 

Now, let us consider an arbitrary canonical transformation 
{q,p) (A,/i) 



(3.23) 



independent of a Casimir coordinate c (not necessarily a point transformation!). The 
advantage of staying inside such a class of transformations is that the clear structure 
of a pencil is preserved and the Marsden-Ratiu projection of the Poisson pencil is still 
trivial. Of course, the most general case of multi-Hamiltonian separability theory takes 
place when one goes beyond the set of canonical coordinates, i.e. when one tries to find 
DN coordinates starting from the pencil written in a non-canonical representation. But 
then the simple structure of degenerated Poisson pencil is lost and the nontrivial problem 
of the Marsden-Ratiu projection for such pencil appears (see for example ref. pl|| ). 

Applying the transformation ( ^.23 ) to Hamiltonian functions (3^) and Poisson matrices 
(|3.12| ) one finds that 



hkiq,P,c) = hkiq,p) +cbk{q,p), k = 1, 



(3.24) 



and 



TTl 




/ 
-/ 



0, 



D{q,p) A{q,p) 
-A^iq,p) B{q,p) 



(3.25) 



where A, B and D are n x n matrices. The nondegenerate Poisson pencil 9x on S gives 
rise to the related Nijenhuis tensor N and its adjoint A^* in (g, p) coordinates in the form 

B )^ ^ °^' = [ D A 



N 



(3.26) 



Obviously, in a real situation we start from a given bi-Hamiltonian chain (|]2|)-(|]2|) 
in canonical coordinates {q,p,c), derived by some method (see for example |13(|-[|25|), and 
trying to find the DN coordinates which diagonalize the appropriate Nijenhuis tensor and 
are separated coordinates for the considered system. So now we pass to a systematic 



derivation of the inverse of transformation (3.23). 

The important intermediate step of the construction of DN coordinates are the so called 
Hankel-Frobenious (HF) non-canonical coordinates {u, v) ||3^ , ||4^ , related to the DN ones 
through the following transformation 



Ui = Pj(Ai, A„), 

n 

l^i = ^UkK-\ i = l,...,n. 

k=l 

In (u, v) coordinates one finds 



(3.27) 



U 
-U^ 



u 



/ 
••• 

V 1 ui 



1 \ 

Ul 
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N 



F 
F 



F 





-FU'^ 



-Ui 


1 







-U2 





1 ••• 





-Un-1 







1 


-Un 










FU \ 

r 


N* 







01 = N 000 = 

\ —r u J 

Moreover the differentials dui, dvi satisfy the following recursion relations [EH], [41 1 



N* o dui 
N* o du2 

N* o dUn-l 
N* o dUn 

N* o dvi 

N* o dV2 
N* o dVn-1 

N* o dv„ 



du2 — uidul, 
dus — U2dui, 

dUn - Un-ldUl, 
—Undui, 

dv2 — uidvi, 
dvs - U2dvi, 

dVn - Un-ldvi, 
—Undvi. 



(3.28) 



Note that vector fields ° dui and 9q o dvi are quasi-bi-Hamiltonian. 

Now we relate the canonical coordinates {q,p,c) to the DN separated coordinates 
(A, /i, c). From the minimal polynomial of N ( |3.26| ) we get 

Uk = Ckil^p), k = l, n. (3.29) 
Conjugate coordinates Vk = Vk{q,p),k = 2,...,n are found from the recursion formula 



(3.30) 



( 3.2§| ) while vi = vi(q,p) coordinate from relations 

{uj,vi}ei^ = 6j^n-k, j = l,...,n, k = 0,...n-l. 
Hence we get 

Vk = i^k{q,p), k = l,...,n. (3.31) 

Eliminating {u, v) coordinates from ( ^) , ( ^2^ and (|^ we derive the desired relations 

Xi(9,P; A,/u) = 0, i = l,...,2n. (3.32) 

Now let us concentrate on a special but important case of point transformation between 
{q,p,c) and (A,^,c) variables. Then 



01 
N 



Aiq) 

-A^{q) B{q,p) 

Aiq) 

Biq,p) A^iq) 



(3.33) 
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where matrix elements of B are at most linear in p coordinates, so coefficients of minimal 
polynomial of N are equal to coefficients of the characteristic polynomial of A. Hence, we 
find the first part of the canonical transformation in the form 



PiW = Vi{Q), i = l,-.,n =^ qi, = iPf^{X), k = l,...,n. (3.34) 
The complementary part of the transformation we get from the generating function 



i=l 



Then, 



dG 



1, ...,n 



Pk = <y5fc(-^>M), k = l,...,n. 



(3.35) 



(3.36) 



At the end of this subsection we introduce the notion of an inverse bi-Hamiltonian 



separable chains. In ref. |2£] it was demonstrated that for each separable bi-Hamiltonian 
chain ( |3.3[ ), ( |3.24| ), ( p. 25 ), ( |3.26| ) in canonical coordinates {q,p,c) there exists a related 
inverse bi-Hamiltonian separable chain 

ttq o dh'^^^^ = 



TTo o dh'^ 



K'^= vT-i o dh[.^i 



(3.37) 



ttq o dh'^ = K[ = 7r_i o dh'2 
= 7r_i o dh[ 

where /i^+i = /iq = c. 



0-1 K„ 

— iT 

-K„ 



= Oq ° dh'^ and 

KiQ^P^c) = hr{q,p) +cb'j,{q), h'^{q) = r = l,...,n. 



(3.38) 



(3.39) 



Notice that in both chains the respective Hamiltonians ( 3.24p , ( p. 39 ) and related vector 
fields differ only by the c— dependent parts. In the case of a point transformation to the 
DN coordinates, when 9i takes the form ( |3.33| ), we get 



{A-^foDoA-^ -{A- 
A-^ 



(3.40) 



and 



o 67^ o 6*0 



A- 



(3.41) 
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Notice that both chains ( |3.3D and ( |3.37| ) given in natural coordinates can be transformed 
to the Nijenhuis bi-Hamiltonian form .61 ) -( ^112 ), where the canonical transformation can 



be derived from the relation Pr{^) = bj.{q), r = l,...,n in the first case and from the 
relations /Oj,(A) = r = 1, ...,n in the second case. 

Consider the set of extended functions 

hr{q,p; c, c) = hr{q,p) + cbr{q) + cb'{q). (3.42) 

They can be simultaneously put into the bi-Hamiltonian and the inverse bi-Hamiltonian 
hierarchies. In the first case c is treated as the Casimir variable and c' as the parameter 
and in the second case c is treated as the parameter and c' as the Casimir variable. 

3.3 Examples 

We shall illustrate the theory presented by a few representative examples. More examples 



can be found in refs. [26|-||2^. In all examples from this section canonical transformations 
between natural and separated coordinates are point like. An example of a nonpoint 
transformation will be given at the end of this paper. 

Example 1. The one- Casimir extension of the Henon-Heiles system. 

Let us consider the integrable case of the Henon-Heiles system generated by the Hamilto- 

nian H = ^pf + + 9i + ^91 ^i- one-Casimir extension reads [22| 



{qi)tt = -3gf - + c, {q2)u = -qm (3.43) 
and belongs to the bi-Hamiltonian chain 



(3.44) 





o d/iQ 


= 








o dhi 


= Ki 


= TTl 


o dh^ 




o dh2 


= K2 


= TTl 


o dhi 









= TTl 


o dh2 



where 



ho = c, 

hi = hi{q,p) + cbi{q) = -pj + -pi + qf + -qiql - cqi, 

2 2 2 (345) 

1 1111 

h2 = h2{q,p) + cb2{q) = -q2PiP2 - -qipl + — g| + -7qfq2 - tCqI, 

2 2 lb 4 4 



/ 








qi 


\q2 


Pi 


\ 










\q2 





P2 






-qi 


-^2 





\P2 


-hl,qi 






-\q2 





-\P2 





~hl,q2 




V 


-Pi 


-P2 


hl,gi 


hl,q2 





/ 



(3.46) 
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The construction of the related inverse bi-Hamiltonian chain, according to the results from 
the previous subsection, gives the following results: 



ttq o (i/ig = 

ttq o = K2 = vr_i o dh'^ 
ttq o dh'i = K[ = 7r_i o 

= 7r_i o dh'i, 



(3.47) 



where 



4 

K = hi{q,p) - — c, 

W h t (3-48) 
1^2 = h2{q,p) + —c, ^ ^ 

Q2 



h'n = c, 



vr-i 












2/92 


I92P2 


\ 








2/(72 


-4gi/gi 


^92^1 - Q'lP2 







-2/g2 





2p2/(?i 






-2/(?2 


4gi/gi 


-2p2/gi 











-^92Pi + giP2 


^2,,, 


^2,,. 





/ 



(3.49) 



and Newton equations related with the natural Hamiltonian h'^ are: 

3 



1 



(9i)« = -3gi - -g2; (g2)tt 



-gig2 



(I2) 



(3.50) 



being just the second well known one-Casimir extension |16] of the Henon-Heiles system 
considered. Both systems ( p. 43 ) and (3.50) can be transformed to the Nijenhuis chain 
j.6|)-( |3lT^) through the respective transformations 



gi = Ai + A2, v\ 



+ 



52 = 2v^-AiA2, P2 
1 



^2\^2 

A2 — Ai ' 
^1 



+ 



1^2 



Ao A 



Ai 



(3.51) 



and 



gi 



g2 



1 ^ 1 
Ai A2 ' 



A1A2 



AiA^i j_ A2A*2 



Ai — A2 



\/— A1A2 



P2 
4, .2 



'^-^(aT 



A?^i 



A2 — Ai 

+ 



A2^2 



A2 A2 



(3.52) 



i^Kf4 + A 



-3 
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Example 2. One-Casimir extension of the Kepler problem in the plane. 

Let us consider the classical problem of a particle in the plane under the influence of the 

Kepler potential and an additional homogeneous field force. The Hamiltonian function 

reads 



hi{q,p,c) 



1 2 , 1 2 Q 



cq2, a = const. 



(3.53) 



cqf, 



(3.54) 



There is a second independent integral of the motion 

, . X 1 2 1 1 ags 1 

h2{q,P,c) = --q2Pi + ^qiPiP2 + -—^=== - - 
2 2 2 + g2 4 

which together with h^ = c allows us to construct a bi-Hamiltonian chain (3.44) with the 
second Poisson structure in the form 

/ PI \ 

P2 

-hi^qi ■ (3.55) 

"^1,92 





TTl 



V 












-q2 

-P2 





\qi 


\PI 

''1,91 



92 




'-1,92 



-Pl -P2 hi 

The inverse bi-Hamiltonian chain ( |3.47| ) is given for functions 

a 4 

h2[q,p,c) 



h[{q,p,c) 



1 2^ 1 2 
1 



V^l + 92 



2 ± I aq2 



+ ^c, 
9f 



(3.56) 



/i3 = c, 

and the second Poisson tensor in the form 



7r_i 



/ 








-492/9? 


2/91 


^92Pi - 91P2 \ 










2/91 





\q2P2 




492/91 


-2/91 





-2pi/9f 


-h' 




-2/91 





2pi/9i 





"^2,52 


V 


-\q2p1 + 91P2 


-\q2p2 


^2,9, 


^2,52 


/ 



(3.57) 



The transformations to DN coordinates for both chains are the following 

A^i , IH 



91 



2v^AiA2, Pl = sj A1A2 



§2 = Ai + A2, P2 



A2 

AlAii A2/i2 



+ 



A2 — Ai 



Ai 

1x2 1 

-Aj/ij + -a 



+ 



A2 A 



Ai 



(3.58) 



and 



91 



92 



a/— A1A2 

Ai + A2 



A?/ii 



+ 



^2/^2 



AiA 



P2 



A1A2 



1^2 



Ai^i 



+ 



A2 A2 



Ai — A2 A2 — A 



(3.59) 
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Example 3. One-Casimir extension of elliptic separable potentials. 



In refs. pq| ,|2S] it was proved that every natural Hamiltonian system 

1 " 1 
H=-Y,pl + V{q) + -ciq,q), 

i=l 



(3.60) 



where (., .) means the scalar product, which admits in extended phase space M 3 {q,p, c) 
the bi-Hamiltonian formulation 



q 

P 
c 



ttq o dhi = VTi o dho, 



(3.61) 



where ttq is a canonical Poisson structure, 



VTi 



/ /ii,p 

-A+^q(g)q \p®q-\q®p 



(3.62) 



A = diag{ai, aj\f), Oj— different positive constants, Hq = c, hi = H + cpi{a), is sep- 
arable in generalized elliptic coordinates. This bi-Hamiltonian formulation generates the 



chain (1^) [25| of commuting bi-Hamiltonian vector fields, where 
hr{q,p,c;a) = hr{q,p;a) +cbr{q]a), 

briq;a) = Pr{a) + - '^Pkia){q,A'' ^ ^q), r = l,...,n-l. 



k=l 



bn(a) = Pn{a)[l - -{q, A ^q)], 



(3.63) 



hr{q,p;a) = ^pfc(a)/i, 

fc=0 

liPj - QjPi 



1 
2 



i=l 
n 



^ — ^ rv ■ — rv ■ ^ — ^ 



i=l 



{H, Ri}no = Oi ^ = •••) and Pria) are Viete polynomials of a. 

On the other hand, according to our procedure, the inverse bi-Hamiltonian chain ( f3.47 ) 
for one-Casimir extension of potentials separable in elliptic coordinates reads 



h'^{q,p,c;a) = hr{q,p;a) + c 
/ 



br-i{q;a) 

bniq;a) 



'n,p 

B-^ -B-^{\p®q-\q®p)B-^ -hn,g 
\ - {hn,pf {K,pf 



(3.64) 
(3.65) 
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where 



B = A-^q^q, 



B 

\B 

da = diag 



9a \B\ + -{daq®daq) \A\ 



d_ _d_ 

' dar, 



dai 
1 



(3.66) 



\B\ = \A\--{q,da\A\q). 

Notice that now, the natural Hamiltonian in the inverse hierarchy is the last one of the 
form 



H' = h[{q,p,c) = ]-j2p' + ^(l^ + 



i=l 



Pn{a)[l - {q,A~\)y 



(3.67) 



The basic example here is the Garnier system with the potential V{q) = ^{q, q)'^ — -^{q, Aq). 
This potential is a member of an infinite family of permutationally symmetric potentials 
separable in generalized elliptic coordinates |17]. The point transformation to DN sepa- 
rated coordinates can be constructed from the relations 

1 '"^ 

p,{X) = p,{a) + -Y,Pkia)iq,^'"^~^'l), r = l,...,n. (3.68) 

fc=i 

But so defined DN coordinates are just the generalized elliptic coordinates Ai, A„ defined 
by the relation 

The proof is given in refs. [^6| and [^. 

More examples of separated systems by the method presented the reader can find in 
refs. [||]-||,l3l-[i|. 



4 Multi-Casimir unsplit chains 

In the previous section a separability theory of one-Casimir bi-Hamiltonian chains was 
reviewed. Here we pass to the generalization of the theory and include multi-Casimir 
cases. This procedure considerably extends the class of separable systems and in general 
covers a new class of chains, i.e. the so-called split chains. In the following section we 
consider the simplest generalization of unsplit multi-Casimir chains related to the extension 
of the separation curve (3.4) of the form 

f{X,p) = hx, ^A = c„A2"-i + ... + ciA" + /iiA"-i + ... + /i„. (4.1) 

Choosing other admissible forms of the separation curve with more then one Casimir, 
one can construct split bi(multi)-IIamiltonian chain, i.e. the chain which splits onto a few 
bi(multi)-Hamiltonian sub-chains, each starting and terminating with some Casimir of the 
appropriate Poisson structure. The work on split cases is still in progress but some results 
are presented in the next section. 
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4.1 Multi-Hamiltonian Darboux-Nijenhuis chains 

In ref. pO| ] the multi-Hamiltonian chain was constructed for the separation curve in the 
form ([4.1|). Actually, the Hamiltonian functions hk found from the system 



f{Xi,fii) = Cn}?^'^ + ... + ciA^ + /iiA^-^ + ... + /in, i = 1, n, 
take the following form 



hr{X,fi,c) = - ^ 



1=1 



dXi Ai 



(4.2) 



(4.3) 



where /9^^(A) = /0^(A) and ^,m = 2,...,n, are defined by the recursive formula 



Notice that Hamiltonians (14. 3| ) are just Hamiltonians (|3.6| ) supplemented with extra terms, 
linear with respect to additional Casimirs Cj, i = 2, ...,n. On the extended phase space 
M 3 {Xi, Xn, fJ-i, ci,...,c„,) functions ([4. 3D and /ii_r-(c) = Cr, r = form 



^ J bi-Hamiltonian chains, each generated by a Poisson pencil vr^^fe-i = vTfc — A'^ Vj 
of order 1 < (A; — i) < n, 

VTj o dh-i = 

TTj o dh-i+i = Ki = TTfc o dh^k+i 

TTi O dh^ij^2 = K2 = VTfe O dh^k+2 

■Ki o dh^i+j = Kj = TTfc o dh_k+j 



TT^k-i O d/l;^ 



< i < < n, (4.5) 



VTj O dh-i^n = Kn = TTk ° dh^k+n 

= TTjfc o dh_k+n+l 

with respect to (n + 1) compatible 3n x 3n Poisson structures of rank 2n 



vr2 



/ ^0 


V 



/ 02 
— T 





V 







\ 



TTi 



Ko K, 







V 
\ 



(4.6) 
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TTr, 



— T 
— T 

-K 



n—l 



where 



A™ 
-A™ 



(4.7) 



Km = {hm,fi, —hm,x)'^ and each Poisson structure tt^ has n Casimir functions: 
Cm+i, Cm+2, ■ ■ ■ ,Cn,hn,--- ,/in-m+i- Moreover all functions hr{X,fJ,,c), r = l,...,n are in 
involution with respect to an arbitrary Poisson tensor vr^, k = 1, ...,n. 



Now, we integrate equations of motion from the hierarchy (|4.5| ) solving the Hamilton- 
Jacobi equation for Hamiltonians ([4.3|) 



hr{X. 



dW 



pti(A)/fe(Afc,9VF/aAfc) 



Or- 



(4.8) 



k=l 



i=l 



First we demonstrate the separability of this equation. Taking the generating function 
W{X, a) in the form W{X, a) = XlILi ^ii^i^ and the following representation of Pf, ^ 



n+fc-l 



1=1 

eq.([4.8D turns into the form 



i=l 



I 

•r-1 



n+k—1 



(4.9) 



n k 



k=l 



p':-i{ m{^k,dWk/dXk)-U=i c.A^ 

At 



n—l+i] 



Or- 



Applying relation ( ^.141) we get the solution of eg. ( 4. 10 ) in the form 
f{Xk,dWk/dXk) = g{Xk), k=l,...,n, 

where 

9{C) = cne^"-^ + ... + ciC + ai^"' + ... + an-ii + an. 



(4.10) 



(4.11) 



(4.12) 



Hence, W{X,a) can be obtained by solving n decoupled first-order ODEs ( 4.11| ). For 
example, if 



f{Xi,^ii) = ip{Xi)f{ni) + V'(Ai), 
then we obtain 



k=l 



(4.13) 



(4.14) 
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In new canonical variables Oj, bi — , the Haniiltonians /ij.(A,/i,c) become hj- — Oj. with 

S^-eT (/-')' S^e. ("5) 



where (/~^)' means the derivative of/^^. As in the new coordinates each hr generates a 
trivial flow 

i^X = "l^; = 0' (b^)^r =g^= hr, = 0, (4.16) 

hence 

hi = ti + const. (4-17) 

Combining ( |4.15|) with (|4.17[ ) we arrive at implicit solutions for the trajectories Xi{tr), 
with respect to the evolution parameter tr in the form 

Y^J {f-^y^dC = di,rtr + const, i = l,...,n. (4.18) 

4.2 Multi-Hamiltonian chains in arbitrary canonical coordinates 

Let us introduce arbitrary canonical coordinates {q,p,c) related to the Darboux-Nijenhuis 
coordinates (A, /U, c) trough some canonical transformation 

9fc = Cfc(A,Ai), pk = rif,{X,fi), k = l,...,n. (4.19) 

Applying the inverse of this transformation to Hamiltonian functions (4.3) and Poisson 
matrices (4.6) one finds that 

71 

hr{q,p,c) = hr{q,p) +^Q6i,.r(g) (4.20) 
1=1 

and the nondegenerate part 0^ of rank 2n of each 7rm(also implectic) takes now the form 

a Arm „ /I ( Dm {q,p) Am{q,p) \ _ 1 ^ M o1^ 

^"^"^ °^'=[-Aliq,p) Bm{q,p))^ m=l,...,n. (4.21) 

Conversely, if we have a multi-Hamiltonian chain in {q,p, c) coordinates and the Nijen- 
huis tensor 



Niq,p) 



Di{q,p) Ai{q,p) \ ( l \ ^ _ ( Ai{q,p) -Di{q,p) 



-Al{q,p) B,{q,p) J'\-I J y B,{q,p) Ajiq,p) 

(4.22) 



is nondegenerate and has n distinct eigenvalues Aj each of multiplicity 2, then the canonical 
transformation (4.19) transforms a given chain to the one considered in the previous 
subsection. 
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The admissible reductions of the number of Casimir variables are the following. For 
arbitrary 1 < m < n, let Cj 7^ 0, 1 < i < m and q = for m < i < n. The first 
m Poisson structures vTj, i < m survive the projection (A, ^, ci, c„) £ M ^ M 3 
{X, fi,ci, ...,Cm) and we have still ("^2"^) bi-Hamiltonian chains ( |4.5| ). In the limit ci = 
... = Cn = 0, the systems considered lose the bi-Hamiltonian property, turning into the 
quasi-bi-Hamiltonian systems on a symplectic manifold M 3 {q,p), being still separable 
and integrable by quadratures. Moreover, because of the property (^), each of the 
multi-Hamiltonian systems considered on a Poisson manifold M 3 {q,p, c) has a quasi-bi- 
Hamiltonian representation on a symplectic leaf S of ttq (dim S = 2n) fixing the values of 
all Ci. 

4.3 Examples 

The theory extended in this section will be illustrated by several representative examples 
of already known as well as new multi-Hamiltonian systems. 

Example 4. Stationary t2—flow of dispersive water waves. 

The Hamiltonian functions and Poisson structures in Ostrogradsky variables are as follows 
11: 

hi{q,p,c) = -4pip2 + 5q2pl - ^qiq^ - 92 - + ^q2Ci + {^qi + ^gi)c2, 

h2iq,P, c) = qipj + 4:q2PiP2 - ^qjpl - 2pi + ^9i92 " ^9?9| - ^9? + ^^^92 

13 13 

+ (2^1 + g92)ci - {^qiq2 - Y^92)C2, 



vro = 



TTi = 



/ 








1 





\ 





















1 















-1 



























-1 












































V 














J 










/ 
















1„ 15 „2 
"2^1 - -sl2 




\ 
















1 


Q2 









hi 






-1 







-Pl 


-^1,91 







^ 8 


9i 


-q2 


Pi 


















-hi,p2 


h 


1,91 


^1,52 








V 























0/ 



7r2 = 



/ 





3r,2 1„ 

8^2 - 2^1 


-39192 - ifgi 


^2,pi 


^1,P1 


\ 








-h2 


1^ 7^2 
-2^1 - 8^2 


^2,p2 


^1,P2 




-§92 + ki 


h2 





h2Pi 


-/l2,gi 


-^1,91 




hiq2 + ifgi 


hi + h2 


-h2Pi 





-h2,q2 


-^1,92 




— ^2,pi 


— h2,p2 


^2,qi 


^2,52 










V -hi,pi 


-hi,p2 


^1,91 


^1,92 








/ 
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we have three bi-Hamiltonian chains 










ttq o dci = 









dci = 









ttq o dhi = 


Ki 


= TTi dCl 




dhi = 


Ki 


= T^2 


dC2 


TTo O d/l2 = 


K2 


= TTi d/li 




dh2 = 


K2 


= T^2 


dci 







= TTi (i/l2, 









= vr2 


dhi 


VTi O dC2 = 

















TTi O dci = 


Ki 


= 7r2 dC2 












TTi o dhi = 


K2 


= 7r2 dci 

















= 7r2 d/ii. 













(4.23) 



The canonical transformation to the Darboux-Nijenhuis coordinates reads 

qi = -(3A? + 3Ai + 4AiA2), 
92 = -2(Ai + A2), 
1/^2-/^1 



Pi 

P2 



2 Ai - A2 ' 
1 Ai(3^2 - - ^2(3^1 - 2^2) 



2 Ai - A2 

where now /i^, r = 1, 2 take the form ( |4.3| ) with 

/,(A,,;uJ = 2Af-iM?, 



/5i,i = Pi'/3i,2 = />2 O and 

/52,1 = —'^l — -^2 ~ A1A2, 
1^2,2 = '^l'^2(Al + A2). 



(4.24) 



Example 5. Two-Casimir extension of the Henon-Heiles system. 

In Example 1 one-Casimir extension of the Henon-Heiles system was considered in the 
form 

{qi)tt = -3g? - ^^2 + c, {q2)tt = -qiq2, 



with two constants of motion ( 3.45D and the related transformation to Darboux-Nijenhuis 



coordinates ( 3.51j ). For a two-Casimir extension we get immediately 



P = -[(Ai + X2f - A1A2] = -{qf + \ql), 

1 4 (425) 

(32,2 = AiA2(Ai + A2) = -^qiqi- 



Hence 



hi = H = ]-pI + ]-pI + q\ + \qiql - ciqi - {qf + \ql)c2, 

11 11 1 I (^-'^^ 

h2 = ^q2PiP2 - ^qipl + Y^qt + -^qlql - -^qhi - -^qiqc2, 
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where the Newton's equations related to the energy H are 

{qi)tt = -3qj - ]^ql + ci + 2giC2, fe)** = -^192 + ^g2C2- 
This is tri-Hamiltonian system with the following Poisson structures 



(4.27) 



VTi 



T^2 



1 ° 





1 








\ 











1 








-1 




















-1 
































V 














) 








qi 


\q2 























hl,p2 













-qi 


-^2 


\P2 















-^2 


-ip2 





-^1,92 













—hi,pi — 




^1,92 














V 






















/ 





qf + hi 


\(ll(l2 


^2,pi 


hi,pi 


\ 










\qlq2 


4^2 


^2,p2 


^1,P2 






_„2 _ 1 „2 
Ql 4^2 


-^9192 





^9lP2 


-h2,qi 


-hl,qi 






-^9192 


492 


-^9lP2 





-^2,92 


-^1,92 






—h2,pi 






^2,92 










V 




~hl,p2 




^1,92 








/ 



The first two of them come from one-Casimir extension and the last one was constructed 
according to formula (4.21). Notice that again we have three bi-Hamiltonian chains ( [4.23 ). 

Also in Example 1 the inverse one-Casimir extension of the Henon-Heiles system was 
considered in the form 



1 tt 



-3qt 



1 



Q2, 



2 tt 



-qm 



<l2 



with two constants of motion (3.48) and the respective transformation to the DN coordi- 
nates ( tj.52| ). The two-Casimir extension we get by adding new terms 



/3 _[(Ai + A2)2-AiA2] = 4-^> 

qi Q2 

^2,2 = AlA2(Ai + A2) = 

Q2 



to the constants of motion ( |3.4S| ) 



h[{q,p,ci,C2) = h[{q,p,c = ci) + C2, 

Q2 

h2iq,p, ci, C2) = h2{q,p, c = ci) + ( — - ) C2, 

\Q2 Q2 
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where the Newton equations related to the energy h'^ are 



{Qi)tt = -Sqf - ]rql - ^C2, fe)^ 



8 64gi 

-9192 qCi H F-C2. 

% Q2 



Again this is inverse tri-Hamiltonian system 



TTO 


o dci = 













dci = 













o = 




= vr 


~i 


o dci 




ci/i2 = 


K'2 


= vr 


-2 


(iC2 




dh'i = 




= vr 


-1 


o dh'2 




d/i'i = 


K[ 


= vr 


-2 


dci 









= vr 


-1 


dh'i , 









= vr 


-^2 


dh'2 



(4.28) 



vr_i o dc2 = 
vr_i o dci 
vr_i o d/i2 = K'l 




7r_2 o dc2 

71-2 ° C^Cl 



7r-2 



•dh'o 



where the Poisson structure vr_iis given by (3.49), with additional last row and column 
with zeros, while the new third Poisson structure vr_2 reads 



7r-2 











4/9i 


-8g/gi 
(16g2+4gi)M 




h' 


\ 












h' 


h' 








8qi/ql 





-8qiP2/qi 




-h'2,,. 








-{16qj+4ql)/qi 


891^2/92 

















hi,,. 












V 


-h'2,p. 






h2,,. 








/ 



Till now all examples presented were bi(multi)-Hamiltonian Stackel systems [27|, i.e. 
the systems with all Hamiltonian functions quadratic in momenta in DN coordinates: 
/(Aj, /Xj) = ip{\i)fif + ipiXi). Here we present the first example of non-Stackel system. 

Example 6. m-Casimir extension of the relativistic n-body problem. 
Consider the Hamiltonian dynamical system with the Hamiltonian given by 



H 



E 



A,; 



(4.29) 



where ip^ are arbitrary smooth functions and a is an arbitrary constant. The corresponding 
dynamical system takes the form 



(A 



i)tt 



(Ai)t(Afc)t 



Ai — Afc 



1, ...,n. 



(4.30) 



which depends explicitly on velocities. The derivation of formula ( [4.30|) is given in ref. 
|26|. Notice that equations ( ^.3[1| ) do not depend on cp^ functions, hence the dynamics is 
not influenced by ip^ terms. 

Dynamics ( [4.30| ) is a special case of the integrable relativistic n-body problems intro- 
duced by Ruijsenaars and Schneider [43|. Now, comparing ( [4.29| ) with ( |4.3| ) one immedi- 
ately concludes that (A,/i) is a Darboux-Nijenhuis chart for the Hamiltonian H, and as a 
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consequence, system ( 4.3C| ) is quasi-bi-Hamiltonian and separable, with the solution given 
by the implicit formulae 

-TP: dS, = Si,it + const, i = l,...,n, (4.31) 

where g{^) = a„ + a„_i^ + ... + oi^"""^. This fact was noticed for the first time by Morosi 
and Tondo [^]. Notice, that trajectories Xi{t) do not depend on the factors, as 

expected, so without a loss of generality one can put fi{Xi) = 1- 

The system (|4.3C1| ) can be naturally extended to an m-Casimir one with the Hamiltonian 

^ = E^e"''"+E^AiW' l<m<n (4.32) 

1 i 1 

1=1 j = l 

and the related Newton equations of motion 

k^i * " j=i * 
The dynamical system ( |4.33| ) has n constants of motion 



hr{X, ^, c) = - 2^ — — + 2^ Cj/3j, ^(A), r = 1, n, 
1=1 * * j=i 



{m + 1) Poisson structures (|4.6| ) and the solution given by implicit formula ( [4.31| ), where 
now 

g{0 = an + an-ii + ... + aiC~^ + ciC + ■■■ + c^C^"^-^ ■ 

The one-Casimir extension, which is bi-Hamiltonian, has the following Newton equations 
of motion 

l\ \ {\)t{Xk)t , /\ N -1 

(Ai)tt = 2> — Vac{Xi)t, i = l,...,n. 

k^i ^'-^^ 

The two-Casimir extension, which is tri-Hamiltonian, has the Newton equations in the 
form 

{Xi)tt = ^ X] ^'^\.^^\^ ^ aci{Xi)t + ac2 f Aj + ^ Afc j {Xi)t, i = 1, n. 

k^i * \ k=l / 



Example 7. (m + 1)- Hamiltonian formulation for elliptic separable potentials. 
In Example 3 we presented bi-Hamiltonian systems separable in generalized elliptic coordi- 
nates. In the following example we extend the result onto appropriate multi-Hamiltonian 
chains. 
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According to the theory presented in this section, let us generalize the Hamiltonian 



system (3.60) to the form 

^ m 

H{q,p,c) = -{p,p) + V{q) + ^Ckbk,iiq), 



m = 1, n 



(4.34) 



k=l 



being multi-Hamiltonian and separable. The few first bk,iiq) functions are as follows 



1 1 
b2,i{q) = ^{q,Aq) - -(g,g)^ 

hAq) = ^{q,q){q,Aq) - ^{q,A^ 



(4.35) 



1 



iq,qf,... . 



For example, the three Poisson structures of the system ( 4.34| ) with two Casimirs read 



/ / \ 

-10 



V / 

/ 



VTi 



1 hi^p \ 
-A+^q0q \p®q-\q®p 



1T2 









-{A 



(^qY' 



V 



-{h2,r,? 
-{h 



{hl,qf 

/ 

{A-^q(S)qf 

{A-^q^ q){y (g)q-^q(g)p) 
+ (ij<Og- y(g)p){A- \q^q) 



''2,q 










/ 



The functions hr{q,p,c), forming three admissible bi-Hamiltonian chains, are given by 



formulas ( 3.63 ), where now 
1 " 



r-^Z^^i Ki + Cibi^riq) + C2b2,riq)- 
i=l 



5 Multi-Casimir split chains 

In the following section we extend the results from previous sections onto the so-called 



split bi(multi)-Hamiltonian chains |31]. Potentially, the variety of such systems is much 
richer than the one of unsplit chains, but still less recognized. The reason is that systems 
from this family are generally non-physical in the sense that are either non-Stackel, or 
even if are of Stackel type, the underlying Stackel space is never flat (at most conformally 
flat). 
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Let M be a (2n + k) dimensional manifold endowed with a linear Poisson pencil vr^ 
TTi — Attq. We suppose that it admits k polynomial Casimir functions 



3=0 



a = 1, k 



(5.1) 



with n = ni + ... + 71^. From beeing Casimir of a pencil it follows that the set {^o"^}a=i is a 
set of Casimirs of ttq, while {/in"'*}^^^ is a set of Casimirs of tti, respectively. In canonical 
coordinates {q,p,c) /iq"^ = Cq,, a = l,...,k, 



( Oo ... \ 









V 



VTi 



(1) 



,(fc) 



(5.2) 



where -fcj*^ = Oodh^^' and ^Oi^i given by ( 3.25| ). Now, we are looking for a separation 
curve in the form 



(a) 



(5.3) 



where "da are admissible functions of A and ^. Further on we concentrate on particular 
two-Casimir cases, but it will be enough to make some insight into the theory. 
Let us start from a separation curve for the unsplit two-Casimir case 



/(A, ^l) = hx, hx = C2A"+i + ciA" + /iiA"-^ + 



hr, 



(5.4) 



An arbitrary shift of ci Casimir variable along the polynomial hx leads to a separation 
curve 



/(A, fi) = C2A"+^ + /iiA" + ... + /i,A"-^+^ + ciA"-' + /ii+iA"-^-i + ... + hr. 



(5.5) 



of some split bi-Hamiltonian chain. This is the case ( ^.3| ) with k = 2, ni = n — i, n2 = 

hi, I = 1, i and h^^^ = hi^j,j = 1, n — i. We illustrate 
the situation for i = I. The following results were obtained. For a separation curve in the 
form 



i, = 1,^2 = A"+i~\ /if ^ 



/(A, /i) = C2A"+^ + /iiA" + ciA"-i + /i2A"-2 + ... + /i„ 
Hamiltonian functions hr,r = 1, ...,n read 

hr{X,fl,C) = hr{X,fl) +7^^(A)ci +72,r(^)c2, 

where 



h. 



(5.6) 



(5.7) 
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1=1 

n 



\k=l / kj^i 

7i,i(A) = — , 7i,r(A) = — , r = 2,...,n, 

Pi Pi 

/^^ P1-P2 PlPr+l-p2Pr ^ 

72,i(A = , l2,rW = ' r = 2,...,n, 

Pi Pi 

a* (A) = OriXi = 0), ar{\) = P2,rW = Pr+1 " PlPr- 

Two compatible Poisson structures 



VTi 





0/0 


\ 












-/ 





























v 





) 















A 


h2,^l 




\ 




-A 





-h2,x 


-/ll,A 






-ih2,,f 


ih2,xf 










v 


-ihi,,V 


{hi,xr 








/ 



(5.8) 



(5.9) 



give rise to a linear Poisson pencil 7^;^ = vri — Attq, which acting on its Casimir h\ generates 
a bi-Hamiltonian chain 



VTA o dhx 



TTq dC2 = 







ttq dhi = 


Ki 


= VTi ciC2 







= TTi ci/l^ 


TTq dci = 







TTq dh2 = 


K2 


= TTi dCl 


TTo dhn = 


Kn 


= VTi dhn 







= TTi ci/l„ 



(5.10) 



which splits onto two sub-chains, each starting and terminating with a Casimir of an 
appropriate Poisson structure. 

Example 8. The case of three degrees of freedom. 

We construct a system comparable with the Newton representation of the 7th order KdV 
[2|1, [||]. Thus let us take 



1 



f{fii,Xi) = -fi] + 16\J, i = 1,2,3 



(5.11) 
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and the point transformation generated by relations 



gi = Ai + A2 + A3, 

1 ' \2 , \2 , \2 



^(A^ + + A^) + (A1A2 + A1A3 + A2A3), 



(12 

z 

?3 = 2'''^^ — — A3)(A2 — Ai — A3)(A3 — Ai — A2). 
Then, from (^■^), we arrive at the following Hamiltonians in natural coordinates 

,92^3 



h 



1 = ^ i\pl + PlPS ) + 10^193 + 



+ 



2' 

1 q2 



i - Wqfq2 + 3qt - 4- 



1 

— Cl 



2qi 
1 f2q2_ 
Qi 



qi C2, 



3gi ) P2 + 293^3 + 4 



1 f2q2 



qi 



qi ] piPs - -jPiP2 



:q2P2P3 



- "^qmqs - qlql + ^91^2 + 2g| + ^g?g3 + ql 



qi 



1 



1 



qiq2 



3 1 92 

2'^=^+8^i + 2^,'^^' 



1 2 1 



93 \ ^2 , 1 2 2 , 1 1 93 , 



92 P2 + o^2P3 + 791P1P2 - 7— PlP3 

91/ 8 4 4gi 

- ^ (2^3 + 9192) + -9i93 - 9i9293 - 29193 + 2^i92 - -9i92 " 

^ 9293W 93\ 1.2 ^, 2 ,9293. 

+ ^ + ^ (^92 - -J Cl + -{q,q2 - q,qs + 2q, - 2--)c2. 



9192 



(5.12) 



Hamiltonian split chain (^.10|), 


where 


now 
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^l,Pi 
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^1,P2 
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l92 


^91 


^2, Pi 


^1,P3 


-^91 


-I92 


-93 
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^P3 


-h2,qi 


-hl,qi 
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-^2 


-^P2 








-h2,q2 


-^1,92 





1 
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-^91 


-lP3 








-^2,93 


-^1,93 


— ^2,pi 


—h2,p2 


-/l2,p3 


^2,qi 


^2,92 


h2,q:i 








— ^l,pi 


-^l,P2 


— /ll,P3 


^1,91 


hl,q2 


^1,93 









and are conformally related to these of first Newton representation of 7t/i order stationary 
KdV, i.e. in separated coordinates the metrices of respective Stackel spaces are conformally 
related. 

Another admissible form of the separation curve with two Casimirs, leading to bi- 
Hamiltonian split chain of non-Stackel type, reads 



/(A, fi)=fi (c2A"-* + /iiA"~*-i + ... + K-^) + Cl V + K-i-iX'-^ + ... + K. 
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It has the form ( ^.3D with k = 2, ni = i, n2 = n — i^^Di = 1,^2 = A*, = hi,l = 
1, n — i, and h^^^ = hn-i+j, j = 1, Here we concentrate on the case i = n — 1 

fix, fi) = ^(C2A + hi) + ciA"-i + /i2A""2 + ... + hn. (5.13) 

The Poisson pencil and the chain are the same as in the previous split case, i.e. are of the 
form (|5.9|) and (|5.1C| ), but hr{X,fi,c), being the solution of the system 

f{Xi,fMi) = fii{c2Xi + hi) + ciA""^ + /i2A""^ + ... + /in, i = 1, n, 

are of more complicated form and we omit here the explicit formulas. 

We illustrate the case ( |5.13| ) and n = 3 with the example where the transformation to 
DN coordinates will be of non-point nature. 



Example 9. The 4th order Boussinesq stationary flow. 

The Hamiltonians and compatible Poisson structures for the Ath order Boussinesq sta- 
tionary flow, written in generalized Ostrogradsky canonical coordinates, are the following 



hi= - 27p2P3 - 27p2 - 9pl + 2qlp3 + 3qlp2 + q2Pi - -qiql - -9192 - ^9i 
1 1 

+ -^qici + -93C2, 

2 / N 4 4 o 8 o 16 o 

h2= - 9piP3 + qipi + 2qi{q2 - 2^3)^2 + -^qmPs - ^93 + ^9i93 - ^9i92 

+ ^9i + ^929i - ^9^93 + ^(293 - 92)ci + {3p2 - ^9?)c2, 
/13 = - 54pl - 27p2pl - 81^2^3 + 3giPi + 18^1^2 + 291P3 + M2V1P2 + 3(^2 - 93)^51^3 
+ l5qlp2P3 - 9i(93 + -^q2)pi - i^qf + -9i9293)P2 + (29193 - g9i92 

4 16 4, 32 o 32 fj 16 o 3 4 o 2 

- 9919293 - ^9i)P3 + ^919293 + ^qi + ^91 93 + ^9i92 

4,2 I3 2 o 1. 40 lol 

- ^^293 + ^9293 + ^9293 - ^93 + (^9l + g93 - g9293 - 2^1^32 

- 9iP3)ci + (91P1 + 93P2 - ^9i93 - ^9l92)c2, 



vro = 
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TTi 
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— ^1,P3 


^1,91 


^1,92 
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/ 



where ^ = 9p2 + 3p3 - 9i, -B = ^9 i92 + ^9i93 - 5P1, = -|r9i92 - ^9i93 + 5^1- 

/,From the minimal polynomial ( 3.20| ) of related Nijenhuis tensor one finds the first 
part of transformation to the HF coordinates (li, t;) 



1*2 



93 



:92, 



393 + ^9i - 391P2 - 9iP3 



9 



9293, 



(5.14) 



«3 = ^93 + ^9l93 - 9193P2 - ^9193P3 + ^9i92 - ^9iPi - ^^293, 



while the system ( 3.30 ) with the second chain ( |3.2q ) give the second part of the transfor- 
mation 



Vl 
V2 



91 

92 - 293 
91 



(5.15) 



^^3 = 3p3 + 6p2 



I93 , I9293 4 2 



+ 



9i- 



391 3 91 9 

On the other hand the relation between HF and DN coordinates reads 
Ml = -Ai - A2 - A3, 
U2 = A1A2 + A1A3 + A2A3, 
M3 = — A1A2A3, 



(5.16) 



Hi = ViXj + V2\i + V3, 



1,2,3. 



Hence, eliminating the HF coordinates {u,v) from the system ( |5.14D -( ^.16| ) we arrive at 
the explicit relation between DN coordinates (A,/i) and natural coordinates {q,p). Unfor- 
tunately, the formulas are too long to quote them in the text. In DN coordinates the two 
Poisson structures take the form (|5.9D and the related separated curve 



A"^ = ^(c2A + hi) + ciA^ + /i2A + hs. 



(5.17) 



The implicit solutions of the system can be obtained by solving the three decoupled first 
order ODEs 



dWj 



{c2\i + ai) + {\f + ciXf + 02 Aj + 03), i = 1, 2, 3. 



(5.18) 
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6 Concluding remarks 

In this review paper we presented a multi-Hamiltonian separability theory of finite dimen- 
sional dynamical systems, in the frame of the set of canonical coordinates. It reveals the 
important fact that the multi-Hamiltonian property of considered system is closely related 
to its separability. Actually, we presented the constructive method of finding a separation 
coordinates once having a bi(multi)-IIamiltonian representation of the underlying dynam- 
ical system, written down in some natural canonical coordinates. There is still an open 
question: whether each bi-Hamiltonian chain with sufficient number of constants of mo- 
tion guarantees a separability of underlying Hamiltonian systems? Or, in other words: 
whether the arbitrary degenerate Poisson pencil, written in noncanonical coordinates, can 
be restricted to the nondegenerate one on a symplectic leaf of one of the Poisson ten- 
sors from the pencil? The second important problem is how to perform effectively such 
a Marden-Ratiu projection if it is possible? Some progress in this direction was made 



recently |41] but this part of separability theory still requires further investigations. 
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